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Strong  Optical  Scattering  in  Atmospheric  Turbulence 
Richard  A.  Elliott,  J.  Richard  Kerr,  and  Philip  A.  Pincus 


Abstract 

^A  versatile  and  useful  facility  for  simulating  the  effects 
of  atmospheric  turbulence  on  optical  propagation  is  described  and 
the  relevant  system  parameters  characterized . The  scattering 
medium  is  a turbulent  liquid  (ethanol)  with  the  turbulence  created 
by  unstable  convection  generated  by  a strong  vertical  thermal 
gradient.  Measurements  of  the  structure  function  and  the  spatial 
spectrum  of  the  resulting  refractive  index  fluctuations  are  pre- 
sented and  compared  with  the  theoretically  predicted  forms.  The 
effects  of  this  scattering  medium  on  laser  beam  propagation  are 
determined  and  compared  to  the  first  order  Rytov  theory.  In  par- 
ticular probability  density  functions,  moments  and  spatial  co- 
variance  functions  of  the  irradiance  resulting  from  propagation 
through  the  system  with  a variety  of  turbulence  levels  and  path 
lengths  are  presented. 

The  results  of  an  investigation  into  generating  a controlled 
scattering  medium  by  dispersing  uniformly  sized  transparent  spheres 
in  a medium  with  a slightly  different  index  of  refraction  is  also 


reported . 
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I.  Introduction 

A one-year  program  in  modeling  optical  scattering  from  atmos- 
pheric turbulence  through  the  use  of  a controllable,  repeatable 
laboratory  scattering  medium  has  resulted  in  the  establishment  of 
a versatile  facility.  The  goal  of  this  program  has  been  to  investi- 
gate fundamental  aspects  of  transitional  and  strong  scattering  under 
conditions  which  are  not  readily  achieved  over  a real  atmospheric 
path.  These  aspects,  which  include  the  verification  of  recent 
analytical  predictions,  are  important  for  the  confident  use  of  the 
Green's  function  (Huygens-Fresnel)  formulations  for  complex  sources, 
targets,  and  scenarios.  The  current  re-emphasis  of  shorter-wavelength 
(as  opposed  to  infrared)  systems  gives  these  strong-scattering 
questions  direct  relevance  to  systems  applications. 

The  deleterious  effects  that  propagation  through  a randomly 
inhomogeneous  medium  has  on  an  initially  coherent  optical  signal 
can  only  be  described  in  statistical  terms.  The  probability  dis- 
tribution of  the  irradiance  fluctuations,  its  moments,  and  the 
spatial  and  temporal  covariance  functions  are  examples  of  quan- 
tities which  can  be  measured  and  which  are  important  for  designing 
optical  systems  operating  under  such  conditions.  A considerable 
effort  has  been  expended  in  recent  years  in  attempts  to  understand 
the  relationship  between  the  statistical  properties  of  the  ir- 
radiance and  measurable  characteristics  of  the  random  medium. 

These  characteristics  are  again  statistical  properties  such  as  the 
variance  of  the  index  of  refraction  inhomogeneities  and  their 
spatial  spectrum. 


The  results  of  this  effort  have  been  mixed.  The  Rytov  or 
first  order  theory'*'  has  been  highly  successful  for  weak  scatter- 
ing and  relatively  short  propagation  paths,  where  the  normalized 
variance  of  the  irradiance  o^2  <<  1,  but  it  is  entirely  erroneous 

for  longer  paths  where  the  experimental  results  show  a^2  saturating 
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to  a value  of  order  unity  rather  than  increasing  monotonically . 

More  modern  works  based  on  the  extended  Huygen-Fresnel  principle^  ^ 

8 9 

or  approximate  solutions  for  the  fourth  order  coherence  function  ’ 
have  certainly  elucidated  the  problem  but  the  analytical  expressions 
are  complex  and  it  is  difficult  to  obtain  precise  results  without 
extensive  computation. 

There  are  also  problems  involved  in  obtaining  good  experimental 
data.  Most  of  the  experimental  work  has  been  done  in  the  atmosphere 
because  of  the  potential  impact  on  optical  communication  and  target 
illumination  systems.  However,  it  is  difficult  to  realize  a 
uniform,  stable,  well  characterized  propagation  path  of  sufficient 
length  to  generate  strong  scattering.  For  this  reason  some  model 
system  or  alternative  type  of  scatterer  that  would  allow  experiments 
in  the  laboratory  under  precise  control  is  desirable. 

It  has  been  the  goal  of  this  program  to  obtain  definitive  ex- 
perimental results  under  conditions  allowing  detailed  comparison 
with  the  theoretical  predictions  and  to  determine  empirical  re- 
lationships between  the  statistical  parameters  of  the  irradiance 
and  the  random  medium  for  those  situations  which  still  defy  analysis. 
Our  work  has  proceeded  in  two  parallel  efforts.  One  involves  work 
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with  an  improved  and  refined  version  of  recent  efforts  on  scatter- 
ing by  a turbulent  liquid,  and  the  other  uses  randomly  distributed 
uniform  spheres.  We  will  report  on  the  turbulent  liquid  work  first. 

II . Turbulent  Liquid  Medium 

(i)  Background 

In  the  past  year  there  has  been  significant  activity  in 
this  field.  On  the  analytical  side  Hill  has  made  a more  complete 
study  of  the  spatial  spectrum  of  index  of  refraction  fluctuations 
in  turbulent  fluids.  He  has  considered  not  only  fluctuations  due 
to  temperature  variations  in  air  and  liquids  but  also  due  to  humidity 
fluctuations  in  air  and  salinity  fluctuations  in  water. 

From  our  point  of  view  the  most  significant  of  his  findings  is 
that  the  ubiquitous  Kolmogorov  k-5/ 3 spectrum  does  not  adequately 
describe  the  index  of  refraction  spatial  spectrum  but  must  be  mod- 
ified at  high  wave  numbers.  That  is,  the  -5/3  inertial-convective 
subrange  of  the  spectrum  extends  up  to  < *=  (10n)-1,  where  n = (vVO1'1* 
is  the  Kolmogorov  microscale,  v being  the  kinematic  viscosity  and  e 
the  rate  of  dissipation  of  mechanical  energy.  Beyond  this  point 
and  up  to  < * (Pr)1/2  n_1  there  is  a viscous-convective  subrange 
with  a slope  -1.  For  higher  wavenumbers,  in  the  viscous-dissipative 
subrange,  the  spectrum  falls  off  exponentially.  The  quantity  Pr  is 
the  Prandtl  number,  Pr  * v/D,  where  D is  the  diffusivity  of  the 
quantity  causing  the  index  change.  Hill  and  Clifford  have  also 
calculated  the  effect  this  spectral  change  Induces  in  the  statistical 

properties  of  the  irradiance  and  show  that  it  can  often  be  signifi- 

12 

cant  including  cases  of  practical  interest. 


There  have  also  been  recent,  relevant  experiments  performed 

13 

on  laboratory  generated  turbulence  by  Bissonette  and  Gurvich 
14 

et  al.  Both  groups  heated  a liquid  from  below  and  studied  the 
effect  on  a laser  beam  propagating  through  the  convective  turbu- 
lence. The  use  of  a liquid  rather  than  a gas  is  advantageous  be- 
cause the  change  in  index  of  refraction  with  temperature  (dn/dt) 
is  greater  by  several  orders  of  magnitude,  leading  to  much  stronger 
scattering.  The  Prandtl  number  of  liquids  is  also  higher  (Pr  = 

0.7  for  air,  7 for  water  and  16  for  ethanol)  so  the  one  would  ex- 
pect the  "bump"  in  the  spectrum  predicted  by  Hill  and  Clifford  to 
be  larger.  Unfortunately  the  Russian  group  made  no  attempt  to 

measure  the  spatial  spectrum  of  temperature  fluctuations  directly 

i 

but  assumed  a form  for  the  spectrum  and  determined  the  parameters 
involved  by  measuring  the  thermal  gradient  and  flux  through  the 
cell.  Bissonette  measured  the  temporal  spectrum  of  the  tempera- 
ture fluctuations  and  the  structure  function  of  temperature, 

D^(r)  = <|T(_r  + rj)  - T(r^)|2>,  and  reported  no  disagreement  with 
the  Kolmogorov  spectrum.  However,  since  there  was  no  net  flow  of 
the  water  in  his  scattering  tank,  the  inferred  spatial  spectrum 
is  questionable  because  Taylor's  hypothesis  is  not  valid. ^ Also, 
the  probes  with  which  the  temperature  fluctuations  were  recorded 
could  not  resolve  scale  sizes  less  than  2 mm,  while  scale  sizes 
as  small  as  1 mm  may  very  well  be  present. 

(ii)  Experimental  System 

The  scattering  apparatus  is  illustrated  schematically 
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in  Fig.  1.  The  light  source  is  a 50  mW  HeNe  laser  (Spectra  Physics 
Model  125)  followed  by  a beam  expanding  telescope.  A collimated  beam 
1.2  cm  in  diameter  propagates  through  the  scattering  cell,  a 25  x 25 
x 50  cm  tank  filled  with  ethanol,  at  a height  of  14  cm.  For  multiple 
pass  experiments  mirrors  can  be  inserted  to  provide  up  to  5 passes 
of  the  beam  through  independent  regions  of  the  turbulent  liquid. 

Path  lengths  of  0.5,  1.5  and  2.5  m are  thus  available.  The  light 
leaving  the  cell  is  projected  onto  the  detection  system  of  16  photo- 
diodes in  two  linear  arrays  of  eight  each.  The  effective  inter- 
diode spacing  and  aperture  size  is  controlled  by  the  projection 
optics.  In  the  experiments  described  below  the  aperture  size  ranged 
from  53  to  500  pm.  The  minimum  distance  between  diodes  was  varied 
from  92  to  875  pm  with  the  largest  effective  separation  being  6.13  mm. 
The  signals  from  the  diodes  are  recorded  by  a microprocessor-controlled 
16  channel  14  bit  A/D  converter  and  digital  tape  recorder.  The  data 
collection  system  has  a dynamic  range  of  80  dB  and  the  sampling  rate 
can  be  as  great  as  200  Hz. 

Turbulence  was  created  in  the  scattering  cell  by  maintaining 
an  unstable  temperature  gradient  in  the  ethanol.  The  tank  is  heated 
from  below  by  an  aluminum  plate  resting  on  electrical  heaters  rated 
at  1 kW.  The  temperature  gradient  was  maintained  and  the  mean 
temperature  controlled  by  circulating  thermostatically-controlled 
cold  water  through  a copper  tubing  heat  exchanger  resting  on  a flat 
aluminum  plate  1 cm  below  the  liquid  surface.  In  order  to  prevent 
large  convection  cells  from  forming  vertical  plexiglas  dividers 


Figure  1 

Schematic  Illustration  of 

Turbulent  Liquid  Scattering  Tank  System 


3 cm  high  and  10  cm  apart  run  transversely  and  longitudinally 
across  the  bottom  of  the  tank.  Figure  2 is  a photograph  of  the 
scattering  apparatus. 


The  properties  of  the  turbulent  liquid  which  are  important 
to  optical  propagation  are  the  spatial  spectrum  and  the  structure 
function  of  the  refractive  index  fluctuations.  Since  these  fluc- 
tuations are  due  to  temperature  fluctuations  and  the  rate  of  change 
of  the  refractive  index  with  temperature,  dn/dt,  is  a known  property 
measurement  of  the  spatial  spectrum,  <f>^, ( k)  , and  the  structure 
function,  DT(r)  = <|T(ro  + r)  - T(rQ)  | >,  of  the  temperature  is 
sufficient. 

Structure  functions  were  determined  by  means  of  two  2.5  pm  dia. 
1.5  mm  long  platinum  wire  probes  of  nominal  40  0 resistance,  a Contel 
MT-2  microthermal  analyzer,  and  the  digital  recording  apparatus 
described  above.  The  probe  spacing  was  varied  from  2 to  50  mm. 

The  spatial  spectra  were  recorded  with  a single  probe  (TSI  Model 
1276)  consisting  of  a 25  pm  diameter  glass  fiber  plated  with  platinum. 
The  length  of  the  active  region  of  the  probe  is  250  pm  and  it  is 
presumed  that  it  can  resolve  scale  sizes  of  this  order.  The  time 
response  of  the  probe  was  measured  to  be  ~ 10  ms.  A reversible 
drive  lead  screw  assembly  controlled  by  the  microprocessor  was 
used  to  translate  the  probe  through  the  liquid  and  the  temperature 
recorded  digitally.  The  bandwidth  of  this  recording  system 
(~  100  Hz)  exceeds  by  about  an  order  of  magnitude  that  of  the 
temperature  fluctuations  which  were  observed  in  our  experiments. 
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Figure  f 


(iii)  Experimental  Results 


a)  Turbulence  Measurements 

The  turbulent  temperature  fluctuations  observed  with  this 
system  exhibit  virtually  the  same  properties  observed  in  atmos- 
pheric thermal  turbulence.  These  include  the  'spikey 
randomness'  characteristic  of  small  scale  turbulence  variables 
and  both  spatial  and  temporal  intermittency  often  seen  in  the 
atmospheric  variables.  The  only  significant  qualitative  differ- 
ence we  have  noted  is  the  time  scale  of  the  fluctuations.  In  the 
atmosphere  the  bandwidth  of  the  fluctuations  determined  by  the 
small  scale  sizes  of  the  order  of  millimeters  and  wind  velocities 
of  the  order  of  meters/sec,  is  in  the  hundreds  of  hertz  up  to 
several  kilohertz.  In  the  tank,  the  small  scale  sizes  of  the 
fluctuations  are  also  in  the  millimeter  range  but  the  character- 
istic velocities  are  centimeters/sec  or  less,  thus  yielding  signal 
jandwidths  of  the  order  of  1 to  10  hertz.  Figure  2 displays  a 
single  short  trace  of  the  thermal  fluctuations  observed  with  the 
platinum  thin  film  probe.  In  this  example,  the  probe  was  moving 
through  the  liquid  at  1 cm/sec.  The  fluctuations  observed  are  on 
the  order  of  0.1°  and  indicate  scale  sizes  of  significant  energy 
in  the  range  1 mm  - 1 cm. 

Power  spectra  of  the  thermal  fluctuations  were  obtained  by 
fourier  transforming  recordings  of  the  temperature  like  the  one 
displayed  in  Fig.  3.  Spatial  spectra  were  inferred  from  these 
by  dividing  the  frequency  by  the  probe  velocity  (1  cm/sec)  and 
multiplying  by  2ir.  An  example  of  such  a spatial  spectrum  which 
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Temperature  vs  position.  250  pm  x 25  pm  dia.  cylindrical  thin 
film  probe.  Probe  speed  1 cm/ sec. 


resulted  from  data  obtained  from  60  traverses  of  the  probe  is 
shown  in  Fig.  4.  We  note  that  the  energy  is  down  by  four  orders 
of  magnitude  for  k > 100  cm-1,  i.e.,  for  scale  sizes  less  than 


0.6  mm. 

Once  a spectrum  has  been  obtained  by  this  method  and  stored 

in  the  computer  it  is  a trivial  matter  to  fit  it  to  curves  derived 

from  the  various  theoretical  models.  As  mentioned  earlier  Hill"^ 

12 

and  Hill  and  Clifford  have  determined  that  one-dimensional  spectra 
for  high  Reynolds  number  flows  should  have  an  inertial-convective 
range  where  it  behaves  as  k-5/3,  a viscous-convective  range  varying 
as  k-1  and  a viscous-dissipative  range  which  decays  exponentially. 
Multiplying  our  experimentally  determined  spectra  by  k5^3  or  k-1 
and  replotting  should  therefore  reveal  flat  portions  corresponding 
to  the  inertial-convective  or  viscous-convective  ranges  respectively. 
Fig.  5 displays  the  spectrum  of  Fig.  4 multiplied  by  k5/3.  It  is 
evident  that  the  small  wave  number  portion  of  the  data  is  not  flat 
indicating  that  there  is  not  a well  defined  inertial-convective 
range.  The  form  of  the  one-dimensional  spectrum  expected  at 

12 

higher  wave  numbers  in  a system  with  a large  Prandtl  number  is 
4>t(k)  = Ak-1  exp(-Kii)  where  the  scale  is  related  to  the  Kolmo- 
gorov microscale  n by  n = £ iPr 1 ^ 2/5 . 5 . The  best  least  squares  fit 
of  this  curve  to  the  data  used  in  Figs.  4 and  5 is  shown  in  Fig.  6. 
Fig.  7 displays  a strictly  empirical  fit  of  the  same  data  to  <t>T(x) 

= B(k  + 2k3)-1,  and  Fig. 8 shows  the  best  fit  provided  (^(k)  = ck-1 
exp(-an2<2),  the  Batchelor  spectrum  assumed  by  Gurvich  et  al.^ 


• 


Figure  8 Spatial  power  spectrum  of  temperature  flucuations. 

Solid  line  - Batchelor  Spectrum  <$>(<)  = — exp  (-aic2n2) 

D K 

where  r\/a  = .025  cm. 


LL- 


The  empirical  fit  to  B(k  + b2*3)-1  clearly  agrees  best  with  the 
experimental  data,  a point  made  more  dramatically  evident  in  Figs. 

9,  10,  and  11  where  the  data  points  have  been  multiplied  by  k-1 
exp(-icJli),  (<  + 2k3)-1,  and  k-1  exp(-ouc2)  respectively  and  replotted. 
The  points  should  lie  on  a horizontal  line. 

We  may  compare  the  value  of  n determined  by  the  value  of  i-i 
required  to  fit  the  Hill  spectrum  to  the  experimental  spectrum  with 
that  calculated  on  the  basis  of  the  physical  parameters  of  the  fluid 
and  the  amount  of  power  put  into  the  tank.  The  Kolmogorov  microscale 
is  given  by 


-•(£) 


1/4 


where  v is  the  viscosity  and  £ the  energy  dissipation  rate.  The 
14 


dissipation  rate 


g »TqT 


V 


where  g is  the  gravitation  acceleration,  BT  the  coefficient  of 


thermal  expansion,  q^  the  turbulent  heat  flow  related  to  the  heat- 


ing power  applied  to  the  tank,  the  specific  heat  and  p the  den- 


sity. For  an  input  power  of  900  watts,  tank  area  1250  cm2,  and 
the  parameter  values  appropriate  to  ethanol  given  in  Table  I we 
have  £ = 0.09  cm2/sec3.  This  in  turn  yields  a length  scale  of 
n * 0.054  cm.  This  agrees  quite  well  with  the  experimentally 
determined  value  of  n ■ £jPr1/2/5.5  * 0.04  cm.  It  is  noteworthy 
that  the  length  scale  characteristic  of  the  spectrum  depends  on 
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(k)  exp(an 


• < 


F 

i 
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Table  I 

Properties  of  Ethanol 


density 

P 

0.79  gm/cm3 

index  of  refraction 

n 

1.361 

dn/dT 

-4  x 10"4/°C 

viscosity 

V 

1.5  x 10-2  cm2 /sec 

coefficient  of  thermal 
dif fusivity 

D 

~ 9.4  x 10~4 

coefficient  of  thermal 
expansion 

bt 

0.0011/°C 

specific  heat 

C 

P 

.574  cal/gm°C 

Prandtl  No. 

P 

r 

16 

e-l/4_  Thus,  all  other  parameters  and  conditions  being  equal,  the 
position  of  the  spectrum  of  the  thermal  fluctuations  is  very  in- 
sensitive to  the  power  in  the  turbulence;  for  example,  a 10-fold 
increase  in  heat  flow  qT  yields  a 10-fold  increase  in  the  dissi- 
pation rate  e,  but  only  an  increase  of  the  order  of  1.75  in  the 
value  of  Thus  although  the  characteristic  length  scales  n 

and  £.j  vary  with  the  power  input,  the  variation  is  weak  enough 
that  a single  value  of  is  applicable  to  a substantial  range 
of  heater  power  levels. 

13 

At  this  time  it  may  be  noted  that  Bissonette  reported  ob- 
servations indicating  that  a -5/3  slope  inertial-convective  range 


H 
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spectrum  did  obtain  in  his  experiments  in  turbulent  water.  How- 
ever his  spectra  were  determined  with  a stationary  probe.  We  show 
in  Fig.  12  an  example  of  a power  spectrum  derived  from  stationary 
probe  recordings  of  the  thermal  fluctuations  in  our  ethanol  tank. 
This  data  multiplied  by  f5^3  is  replotted  in  Fig.  13  and  clearly 
exhibits  a flat  portion  at  the  low  frequency  end.  We  infer  from 
this  that  the  -5/3  slope  is  an  artifact  of  measuring  the  spectra 
with  a stationary  probe. 

The  spatial  spectrum  of  index  of  refraction  fluctuation  is 
just  the  spectrum  of  temperature  fluctuations  multiplied  by 
(dn/dT)2,  in  this  case  1.6  x 10-7  °C-2.  Although  all  the  infor- 
mation about  the  turbulent  field  that  is  relevant  to  optical  prop- 
agation is  contained  in  this  function  propagation  theories  are 

often  couched  in  terms  of  the  structure  function  which  may  be 

1 

calculated  from  the  3-dimensional  spectrum,  ^(k)  = “ Yn^  — ^ — » 
by 


Dn(r)  = 8tt  j (l  - -X"r ) <|<n(K)  tc2  d<  (1) 

0 

We  use  the  Hill  spectrum,  = A'  K-1exp(-K£])  with 

A' = (dn/dT)2A  and  = 0.55  mm  as  determined  above,  since  it  leads 

to  a simple  analytic  result  and  since  it  has  some  theoretical 

justification.  The  3-D  spectrum  is  in  this  case 
A'  , 

^(k)  “ YU  K~3  (1  + exp(-Kfci)  and 

D (r)  = C 2 £n(l  + (r/J l^)2)  (2) 

n n 
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S(f)*f5/3 
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where  C^2  = 2a’  is  a parameter  describing  the  strength  of  the  in- 
dex of  refraction  fluctuations  and  corresponds  to  the  usual  struc- 
ture parameter  C^2.  Fig.  14  shows  a plot  of  this  structure  function 

together  with  values  of  <|n(r  + r)  - n(r  )|2>  obtained  from  half- 

o o 1 

hour  averages  of  the  square  of  the  temperature  difference.  This 
exhibits  nice  agreement  over  the  central  portion  of  the  curve. 

The  deviation  at  small  separations  can  be  explained  by  the  in- 
ability of  our  1.5  mm  long  probes  to  adequately  resolve  scale  sizes 
of  this  magnitude. 

The  strength  of  the  turbulence  can  be  characterized  by  the 
constant  C^2  and  this  is  readily  measured  by  determining  Dn(r) 
for  some  convenient  value  of  r.  Fig.  15  is  a plot  of  C^2  deter- 
mined in  this  way,  with  r = 1 cm,  as  a function  of  heater  power. 

The  values  of  C^2  range  over  nearly  two  orders  of  magnitude  from 
5 x 10-11  to  2 x i(T9 . 

b)  Optical  Propagation 

In  part  a)  above  we  have  described  the  characteristics  of  the 
turbulence  in  our  simulation  apparatus  and  have  noted  that  the 
spatial  spectrum  of  the  temperature  or  index  of  refraction  fluctu- 
ations is  different  from  that  observed  in  the  atmosphere.  The  op- 
tical path  length  in  our  system  is  also  much  shorter  than  typical 
atmospheric  propagation  paths.  The  scale  sizes  of  the  turbulent 
eddies  is  somewhat  different  from  those  in  the  atmosphere  and  since 
the  Fresnel  zone  size  in  our  system  is  much  smaller  than  would  be 
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Figure  14  Structure  function  of  index  of  refrac 


the  case  for  a field  experiment  (/XL  = 0.48,  0.83  or  1.08  mm  depend- 
ing on  whether  1,  3 or  5 passes  are  made  through  the  tank)  the 
important  eddy  sizes  are  also  significantly  smaller.  We  therefore 
expect  that  some  features  of  the  radiation  field  may  differ  from 
those  observed  over  long  atmospheric  paths.  In  particular  we  an- 
ticipate a smaller  bandwidth  for  the  intensity  fluctuations  due 
to  the  slower  variation  in  index  of  refraction  fluctuations.  Also 
the  scintillation  patch  size  should  be  smaller  and  aperture  aver- 
aging effects  may  be  important  if  the  detector  size  is  not  suf- 
ficiently small. 

Records  of  the  irradiance  recorded  at  a single  detector  for 
the  weakest  and  strongest  path  integrated  turbulence  are  shown  in 
Figs.  16  and  17.  The  values  of  listed  on  the  figures  are 

2 

the  log-irradiance  variance  calculated  in  the  Rytov  approximation, 
i.e. , 


ini 


8it2k2L  J /l — sin  — 'P  (k)k  d* 

o \ k2L  ' n 


with  k being  the  optical  wave  number,  k the^index  of  refraction 


wave  number,  L the  path  length  and  ^ (K)  = — — (1  + <£i)  exp(-K£,) 

n 4tk  ^ 

with  l i = 0.55  mm.  These  traces  are  qualitatively  similar  to  the 
irradiance  fluctuations  observed  in  the  atmosphere  although  in 
these  laboratory  results  the  time  scale  of  the  fluctuations  is 
significantly  slower. 

The  least  well  resolved  of  the  basic  questions  concerning  the 
statistical  properties  of  theoretical  field  is  the  form  of  the 


probability  density  function  of  the  irradiance.  It  is  generally 
accepted  that  the  PDF  of  the  irradiance  is  log  normal  in  weak 
turbulence  and  exponential  in  an  asymptotic  state  of  very  strong 
path  integrated  turbulence.  Figures  18,  19,  and  20  show  the 
PDF  of  the  irradiance  in  progressively  stronger  path  integrated 
turbulence  conditions.  Also  displayed  in  each  figure  are  the 
maximum  likelihood  fits  to  the  log  normal  and  exponential 
functions.  In  each  case  the  exponential  function  is  a straight 
line.  Fig.  18  is  a weak  turbulence  short  path  condition  and  the 
distribution  shows  reasonable  but  imperfect  agreement  with  the 
log  normal  assumption.  Fig.  19  is  a stronger  turbulence  condition. 
Here  the  normalized  variance  of  the  irradiance  is  near  its  peak, 

* 2.86,  and  while  there  is  scatter  in  the  data,  the  log  normal 
distribution  still  yields  a qualitatively  good  fit  to  the  data. 
Finally  Fig.  20  shows  the  distribution  in  the  case  of  the  strongest 
turbulence  conditions  observed  in  the  system.  In  this  case  the 
dat»  scatter  encompasses  both  the  exponential  and  log  normal  dis- 
tributions except  at  very  low  irradiance  values,  and  we  cannot 
visually  make  a distinction.  We  have  not  been  able  to  reach  strong 
enough  turbulence  levels  to  unambiguously  observe  an  exponential 
distribution.  We  note,  however,  that  in  this  final  example  the 
normalized  variance  of  the  irradiance  is  1.09,  very  close  to  the 
value  of  unity  associated  with  the  exponential  distribution. 

Fig.  21  is  a plot  of  the  normalized  irradiance  variance 

O?  « a?/!*  as  a function  of  turbulence  level  described  by  the 
N 1 


Figure  18  Probability  density  function  of  irradiance.  L = 0.5m, 


Rytov  a*  ^ = 0.19.  Experimental  = 0.72,  = 0. 

N 

A - maximum  likelihood  fit  to  exponential  distribution 


B - maximum  lielihood  fit  to  log-normal  distribution 


Figure  20  Probability  density  function  of  irradiance. 

Rytov  = 210.  Experimental  = 1.09, 

ini  I„  in 

N 

A - maximum  likelihood  fit  to  exponential 
B - maximum  likelihood  fit  to  log-normal 


Normalized  variance  of  the  Irradiance  versus  Rytov  log  irradiance  variance 


Spatial  covariance  functions  were  constructed  from  simultaneous 

recordings  of  the  irradiance  at  the  sixteen  detectors  and  calculation 

of  the  correlation  coefficients  between  the  various  relevant  pairs 

of  detectors.  Different  spacings  and  aperture  averaging  effects 

have  also  been  investigated  by  varying  the  magnification  of  the 

output  beam  before  the  detector  array  and  thus  varying  the  effective 

size  of  the  individual  detectors  and  their  spacing. 

A typical  set  of  covariance  curves  observed  at  one  turbulence 

level  for  three  different  path  lengths  is  shown  in  Fig.  22.  C 2 is 

n 

approximately  8.6  x 10-11  and  the  observed  normalized  irradiance 

variances  are  L - 0.5,  o2  - 0.72;  L - 1.5  m,  o2  *=  2.5;  L * 2.5  m, 

N 

o2  ■ 1.20.  Thus  this  example  encompasses  an  unsaturated  condition, 

N 

a nearly  peak  scintillation  level  condition  and  a strong  saturation 
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re  22  Spatial  covariance  of  the  irradiance.  C2  = 8.6  x 10 


condition.  To  compare  these  to  a first  order  theoretical  co- 
variance  function  we  rescale  the  separation  coordinate  and  replot 
the  data  of  Fig.  22  in  Fig.  23.  The  solid  line  is  that  obtained 
for  a plane  wave  in  the  Rytov  approximation, 


Vd>  ■ Kip  8>2‘2l/*  (l  - 

n \ 


0 

Here  again,  (k)  = (C  2/4ttk3)(1  + <£,-.)  exp(-K£.i)»  has  been  used 
n n 

for  the  spectrum.  The  data  for  all  three  path  lengths  deviate 
significantly  from  the  theoretical  curve  with  the  longer  path  data 
differing  by  the  largest  amount.  The  deviation  of  the  short  path 
data  is  unexpected  and  at  present  we  can  offer  no  explanation. 
Further  investigation  of  this  problem  is  clearly  required.  We  do 
note  however  that  the  greater  path  integrated  turbulence  cases  do 
fall  off  more  rapidly  and  approach  a constant  positive  level  which 
is  the  qualitative  behavior  expected.  Another  example  is  shown 
in  Fig.  24.  Here  C^2  is  approximately  2.4  x 10-10.  With  this 
stronger  turbulence  level,  we  note  the  same  type  behavior  and 
deviation  from  the  theoretical  curve. 

We  next  consider  the  effect  of  receiver  aperture  size  on 
the  covariance  function  measurements.  Figs.  25  and  26  show  co- 
variance  functions  obtained  under  the  same  turbulence  conditions 
as  seen  in  Fig.  24.  Each  figure  is  for  a different  total  path 
length  and  the  three  curves  on  each  figure  are  taken  with  dif- 
ferent effective  receiver  aperture  sizes  ranging  from  0.053  mm 
to  0.167  mm  diameter. 
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Figure  25 


Spatial  covariance  function  of^the  irradiance  for  different 
receiver  aperture  diameters.  C2  = 2.4  x 10-10.  Path  length 
0.5  m.  n 


Spatial  covariance  function  of  the  irradiance  for  different  receiver  aperture  diameters 


For  Fig.  25,  the  turbulence  conditions  and  path  length,  0.5 
meters,  yield  unsaturated  scintillation  conditions.  The  scintil- 
lation patch  size  is  large  and  the  aperture  sizes  used  in  the  ex- 
periment were  all  small  compared  to  this  patch  size.  Thus  aper- 
ture averaging  has  negligible  effect  on  the  observed  covariance 
functions.  In  the  case  of  Fig.  26  L * 2.5  m and  the  scintil- 
lation conditions  are  saturated.  The  smallest  aperture  covariance 
function  shows  a rapid  decrease,  indicating  a small  covariance 
patch  size.  With  increasing  receiver  aperture  size  the  covariance 
function  increases  significantly  indicating  averaging  over  several 
scintillation  patches. 

One  of  the  intentions  of  this  work  was  to  achieve  and  investi- 
gate extremely  high  path  integrated  turbulence  conditions  and  de- 
scribe the  asymptotic  limits  of  the  statistical  properties  of  the 
optical  field  in  strong  saturation.  We  have  found  in  practice 
that  we  could  reach  levels  of  (Rytov)  of  the  order  of  200  with 

reasonable  path  lengths.  While  we  did  not  achieve  as  strong  turbu- 
lence conditions  as  we  would  have  liked,  we  have  been  able  to  measure 
a number  of  important  statistical  characteristics  of  the  optical 
radiation  in  strong  saturation.  These  include  probability  distri- 
bution functions,  normalized  variance  and  spatial  covariance  of 
the  irradiance  as  a function  of  the  turbulence  level  and  path  length. 
III.  Spherical  Scatterers 

We  have  also  explored  the  use  of  a system  of  monodisperse 
spherical  scatterers  as  a model  turbulent  medium.  The  original  intent 


was  to  use  100  ym  diameter  glass  spheres  suspended  in  an  index 
matching  fluid  but  early  in  this  program  it  became  evident  that 
the  uniformity  and  control  required  to  simulate  turbulence  (very 
small  refractive  gradients)  by  using  solid  spheres  in  a liquid 
medium  is  uneconomical  and  basically  impractical.  We  have  there- 
fore pursued  an  alternative  approach  which  is  novel  in  optical 
applications  and  which  should  permit  the  achievement  of  the  de- 
sired scattering  medium.  The  resultant  single-scale-size  scat- 
terers  are  advantageous  in  the  simplicity  they  lend  to  calculations 
e.g.  in  the  multiple  scattering  (saturated  scintillations)  regime. 
An  initial  requirement,  however,  is  to  relate  such  a medium  to 
the  weak  scattering  (turbulence)  regime,  as  follows. 

Let  the  number  of  spheres  of  radius  a per  unit  volume  be 
c and  the  volume  fraction  4>  = (4/3)ira3c  « 1.  The  index  of 
refraction  of  the  scatterers  will  be  taken  to  be  1 + n and  that 
of  the  surrounding  medium  to  be  1.  The  mean  index  of  refraction 
is  then 

<n>  =(l  + m)<(i  + l-  <f> 

* 1 + m<J> 

and  the  variance  in  n is 

on2  = <(n  - <n>)2>  = m24>(l  - 2<|>)  *=  m2<f> 

The  correlation  function  of  the  refractive  index  variation 

Yn(?i  “ r2)  - <(n(rx)  - <n>) (n(r^)  - <n>)> 


J 
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may  be  found  by  calculating  the  probability  that  a rod  of  length 


I 


r = |rj  - r2|  has  both  ends  inside  one  of  the  scatterers.  For  di- 
lute distributions  ($  <<  1)  this  is  easily  effected  and  we  obtain 
the  well  known  correlation  function  for  hard  spheres 

rn(r)  - 1 - <3/2>(i)  + <l/2>(£)3;  0 < t < 2. 

= 0 ; r > 2a 

The  spatial  spectrum  of  index  of  refraction  fluctuations  is 

the  3-d  Fourier  transform  of  a 2y  (r) , i.e. 

n n 
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where  j i (x)  is  the  spherical  Bessel  function  of  order  1. 

The  first  order  theory  of  scattering  from  atmospheric  turbu- 
lence assumes  that  each  scattering  event  is  weak  enough  for  the 
Born  approximation  to  hold  and  that  the  wavelength  X <<  the  size 
of  a turbulent  eddy.  For  refractive  scattering  from  a sphere  of 
radius  a and  index  of  refraction  differing  from  its  surroundings 
by  m this  means  that  Arma/X  <<  1 and  a >>  X. 

This  type  of  scattering  is  commonly  called  Rayleigh-Debye  c-i. 
Rayleigh-Gans  scattering.  Provided  these  restrictions  on  a,  X and 
m are  satisfied  the  above  spatial  spectrum  may  be  freely  substituted 
into  known  expressions  for  e.g.  the  variance  or  covariance  of  the 

O 

log-irradiance  of  an  initially  plane  wave,  i.e.. 
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o2  = 8ir2k2L  f ( 1 — sin  ip  (<)<d< 

ini  Jo  \ <2l  k/  V 

and 

00 

C (r)  = 8ir2k2L  f ( 1 — sin  J (<r)  'I'  (<)<d<  • 

1 0 ' k2L  k / o n 

We  find  then  that  if  the  optical  path  length  through  the  scatter- 
2ra2 

ing  region,  L >>  — - — , the  variance  of  the  log-irradiance. 


o2nj  = 12tt2  m2<paL/X2  and  the  covariance 
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1 - A - r2/4a2 
1 + /l  - r2/4a2 


for  0 < r < 2a  and  Bj(r)  = 0 for  r > 2a. 

This  covariance  curve  normalized  to  unity  is  plotted  in  Fig.  27. 

Experimental  verification  of  this  scattering  theory  requires 
a scattering  medium  consisting  of  uniformly  sized  spheres.  The 
radius  of  these  spheres  must  be  much  greater  than  the  wavelength 
of  the  light  (a  >>  X)  and  the  refractive  index  must  differ  from 
the  intervening  medium  by  an  amount  small  enough  to  ensure  that 
the  scattering  is  weak  (4irma/X  « 1).  These  requirements 
are  relatively  stringent.  It  is  difficult  to  control  the  index 
of  refraction  of  materials  to  better  than  10“4  which  then  implies, 
for  wavelengths  in  the  visible  spectrum,  that  a « 400  pm.  The 
optical  quality  of  such  small  glass  spheres  is  generally  poor  and 


also  the  densities  of  index  matching  fluids  are  usually  much  less 
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than  glass  so  that  maintaining  a distribution  of  glass  beads  in 
a scattering  cell  would  be  a problem.  A more  promising  approach 
and  the  one  which  we  have  pursued  is  to  form  an  emulsion  of  one 
liquid  in  another. 

There  is  a large  class  of  polymers  which  are  soluble  in  water 
and  which,  when  two  dilute  solutions  are  mixed  together,  are 
completely  miscible  but  when  more  concentrated  solutions  are  mixed 
from  two  distinct  phases  with  each  phase  containing  both  polymers 
and  water  but  in  different  concentrations.  The  two  phases  have 
different  densities,  viscosities,  indices  of  refraction,  etc. 

These  two-phase  aqueous  polymer  solutions  have  been  used  by  biolo- 
gists for  separating  different  types  of  cells  and  certain  properties 

19 

of  these  systems  have  been  explored.  However,  their  optical  prop- 
erties have  not  been  extensively  studied.  We  have  initiated  a study 
of  one  particular  system  to  determine  the  feasibility  of  using  such 
systems  for  scattering  experiments. 

The  polymers  we  have  initially  chosen  to  study  are  polyvinyl 
alcohol  (PVA)  and  Dextran.  Fig.  28  shows  the  binodal  curve  for 
these  polymers  dissolved  in  water  at  20°C.  Above  the  curve  the 
solution  separates  into  two  phases  while  for  concentrations 


corresponding  to  points  below  the  curve  only  one  phase  is  present. 


In  all  cases  the  phases  are  optically  clear  with  negligible  ab- 
sorption. A 5%  W/W  Dextran  - 1.5%  W/W  PVA  solution,  for  example, 
separates  into  a top  (10%  vol.)  and  bottom  (9C%  vol.)  phase,  the 
density  of  the  bottom  phase  being  ~ 1.05  g/ml  and  that  of  the  top 
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Figure  28. 


Binodal  curve  for  Polyvinyl  Alcohol  and  Dextran 
in  water  at  20°C. 
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through  it  because  of  the  refractive  scattering  from  the  spherical 
droplets  of  the  dispersed  top  phase.  It  takes  several  hours  for 
the  two  phases  to  separate  after  an  emulsion  has  been  formed. 

The  emulsion  has  been  examined  microscopically  and  it  is  found 

to  consist  of  a dispersion  of  spheres  ranging  in  size  from  ~ 10  pm 

to  300  pm.  We  are  currently  examining  methods  of  producing  a 

monodisperse  emulsion.  One  which  looks  promising  and  will  be 

pursued  involves  subjecting  the  mixture  to  a constant  shear  for 

a period  ranging  up  to  an  hour.  This  produces  droplets  of  the 

20 

dispersed  phase  of  uniform  volume.  The  volume  depends  on  the 
rate  of  shear,  viscosity  and  inter facial  tension.  By  proper 
choice  of  shear  rate  we  expect  to  be  able  to  form  uniform  droplets 
which  can  then  be  used  for  scattering  experiments. 

IV.  Conclusion 

A versatile  and  useful  facility  for  simulating  the  effects 
of  atmospheric  turbulence  on  optical  propagation  has  been  con- 
structed, and  its  parameters  characterized.  The  strength  of  the 
turbulence  generated  by  unstable  convection  has  been  determined 
as  a function  of  power  input  and  the  structure  function  and 
spatial  spectrum  of  the  resulting  index  of  refraction  fluctuations 
recorded  and  compared  with  the  theoretically  predicted  form. 

The  effects  of  this  scattering  medium  on  optical  propagation 
has  been  studied  in  some  detail.  The  probability  density  function 
(PDF)  of  the  irradiance  has  been  determined  for  a variety  of 
turbulence  levels  and  path  lengths  with  the  results  conforming  in 
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a qualitative  way  to  the  theory.  The  PDF  does  not  appear  to  be 
exactly  lognormal  even  at  weak  turbulence  levels  and  short  path 
lengths.  It  does  appear  to  become  progressively  closer  to  an 
exponential  distribution  with  increasing  path  integrated  tur- 
bulence as  predicted  but  in  no  case  is  the  functional  form  un- 
ambiguous. Measurement  of  the  normalized  variance  confirms 
this  behavior  as  it  approaches  unity  as  the  path  integrated 
turbulence  increases  (o  2 = 1.09  when  the  Rytov  logvariance 
o2nI  * 210).  Further  work  in  this  area  should  include  experi- 
ments at  longer  path  lengths  and  also  the  measured  PDF  should 

be  compared  with  the  exponent ial-Bessel  distribution  recently 

21 

proposed  by  Bissonnette  and  Wizinowich. 

The  measured  spatial  covariance  functions  of  the  irradiance 
shows  considerable  deviation  from  the  first  order  (Rytov)  theory 
even  at  weak  turbulence  levels  which  together  with  the  deviation 
of  the  PDF  from  lognormal  may  indicate  that  the  first  order 
approximation  is  invalid  sooner  than  expected.  This  should 
certainly  be  investigated  in  more  depth  and  the  data  compared 
with  predictions  derived  from  the  accurate  multiple  scattering 
theory. Apart  from  this  possible  explanation  there  appears 
to  be  no  reason  for  the  discrepancy.  Aperture  averaging  affects 
can  be  ruled  out  as  they  affect  the  covariance  in  the  opposite 
manner  namely  the  measured  values  would  be  greater  than  the 
true  values.  We  have  also  determined  the  effects  of  the  detector 
aperture  size  directly  (see  Figs.  25  and  26)  and  found  that  the 
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scintillation  patch  size  was  always  greater  than  the  smallest 
effective  aperture  used  (53  pm) • 

There  are  of  course  more  experiments  which  can  be  done  to 
extend  our  understanding  of  the  basic  properties  of  the  scattering 
medium  and  the  effects  on  optical  propagation  particularly  in  the 
multiple  scatter  regime  and  these  should  be  pursued.  However  the 
major  use  for  a facility  such  as  the  one  described  here  will 
probably  be  in  the  investigation  of  such  systems  and  phenomena 
as  speckle  propagation  through  turbulence;  beam  wander-tracking 
and  cancellation  systems;  turbulence  effects  on  radiation  re- 
flected from  glints  and  diffuse  targets;  propagation  of  partially 
coherent  light;  and  pulse  stretching.  Investigations  of  this 
type  will  be  much  easier  to  perform,  and  the  results  will  be 
less  ambiguous  because  of  the  completely  characterized  scattering 
med ium . 
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